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ABSTRACT 

We present optical and X-ray properties for the first confirmed galaxy cluster sample selected by 
the Sunyaev-Zel'dovich Effect from 148 GHz maps over 455 square degrees of sky made with the 
Atacama Cosmology Telescope. These maps, coupled with multi-band imaging on 4-meter-class optical 
telescopes, have yielded a sample of 23 galaxy clusters with redshifts between 0.118 and 1.066. Of these 
23 clusters, 10 are newly discovered. The selection of this sample is approximately mass limited and 
essentially independent of redshift. We provide optical positions, images, redshifts and X-ray fluxes 
and luminosities for the full sample, and X-ray temperatures of an important subset. The mass limit 
of the full sample is around 8.0 x 10""^^ Mq, with a number distribution that peaks around a redshift 
of 0.4. For the 10 highest significance SZE-selected cluster candidates, all of which are optically 
confirmed, the mass threshold is 1 x 10^^ and the redshift range is 0.167 to 1.066. Archival 
observations from Chandra^ XMM-Newton, and 05747" provide X-ray luminosities and temperatures 
that are broadly consistent with this mass threshold. Our optical follow-up procedure also allowed 
us to assess the purity of the ACT cluster sample. Eighty (one hundred) percent of the 148 GHz 
candidates with signal-to- noise ratios greater than 5.1 (5.7) are confirmed as massive clusters. The 
reported sample represents one of the largest SZE-selected sample of massive clusters over all redshifts 
within a cosmologically-significant survey volume, which will enable cosmological studies as well as 
future studies on the evolution, morphology, and stellar populations in the most massive clusters in 
the Universe. 

Subject headings: cosmic background radiation — cosmology: observations — galaxies: distances and 
redshifts — galaxies: clusters: general — large-scale structure of universe 
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1. INTRODUCTION 

Clusters of galaxies mark the largest virialized struc- 
tures in the Universe, and their formation and evolution 
rates depend on cosmological parameters and the kine- 
matics of the dark matter. The cluster mass function 
constrains the normalization of the matter power spec- 
trum, typically expressed as ag, the root-mean-square 
(rms) mass fluctuation on a scale of 8 h~^Mpc. As a 
result, the number of clusters as a function of redshift 
provides a powerful constraint on both the expansion 
history of the Universe and the gravitational growth of 
structure within it (see Carlstrom et al.||2QQ2 , and refer- 
ences therein). Both functions are strongly influenced by 
the matter content Qm and the equation of state of the 
so-called "dark energy" that comprises most of the effec- 
tive energy density in the Universe. Galaxy clusters also 
harbor a significant fraction of the visible baryons in the 
Universe, in the form of a hot intracluster medium that 
leaves an imprint on the Cosmic Microwave Background 
(CMB) radia tion though the Sunya ev Zel'dovich Effect 
(SZE) (S unyaev fc Zel'dovich|ri972[ ). At frequencies be- 
low the SZ null frequency around 218 GHz, clusters ap- 
pear as small (< ImK) decrements in the temperature 
of the CMB. 

Selecting clusters through the strength of their SZE 
signal has key advantages over other methods of selec- 
tion: (1) it is relatively independent of redshift, (2) ac- 
cording to simulations, the SZE flux {Y = J ydCl, where 
y is the usual Compton ^-parameter) should b e cleanly 
rel ated to the cluster mass with low scatter ( |Motl et 
[g^ [2QQ5l [NiiiI|[2QQ6l [Rdd fc Spergel 2006), and (3) ac- 
cording to sfmulations the SZE signal is more robust to 
baryonic astrophysics (cooling, AGN feedback). Thus 
an SZE-selected sample could offer a reliable correlation 
with cluster mass, which is the key cosmological prop- 
erty, particularly at high redshift. However, optical and 
X-ray observations are crucial to confirm the SZE de- 
tections and calibrate the F-mass relation predicted by 
simulations. 

The SZE has been the focus of a considerable num- 
ber of observational efforts over the past three decades. 
So me of the earliest c onvincing measurements were made 
by |Birkinshaw et al.||l984) observing at 20 GHz on the 
40-m telescope of the Owens Valley Radio Observatory. 
Each individual cluster detection required many hours 
of telescope time, typically spread out over several years. 
By the late 1990 's there were still only a dozen or so 
clusters for which reliable detections had been published 



(see the review by Birkins haw|p~999 ) . 

Large-area cosmological SZE surveys have only now 
become a reality, thanks to the development of mm-band 
cameras composed of arrays of order a thousand bolome- 
ters located at the focus of large-diameter telescopes (6 to 
10 m) offering diffraction-limited imaging at arcminute 
resolution (a factor of five finer than the Planck satel- 
lite). Two such experime nts, the Atacama Cosmology 
Telescope (ACT, iFowler et al.,. 2007.) and the South Pole 
Telescope (SPT, |Carlstrom et al. 2009) have surveyed 



hundreds of square degrees of the southern sky at ar 
cminute resolutio ns and are reporting the first cluster 



detections ( Staniszew ski et al. [2009 Hincks et al. 2009 
[Vanderlinde et al. 2010| ). ACT is conducting a wide- 
area mm-band survey at arcminute resolution with :^30 
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Figure 1. The predicted number of clusters more massive than 
4, 6 and 8 x 10^^ Mq (solid lines) as a function of redshift and 
normalized to unity from the mass function of Tinker et al. (2008}. 
For reference we show the 50% and 90% fractional completeness 
levels as dashed lines. The inset shows the normalized differential 
dN/dz distribution over the same redshift range and mass limits. 

/iK s^/^ array sensitivity to efficiently search for massive 
galaxy clusters using the SZE. The first ACT SZE clus- 
te rs de tected in the 2008 observation season are reported 
in iHincks et al. (2009) and the first analysis of the ACT 
148 GHz CMB pow er spectrum reachin g to arcminute 



scales is presented in Fowler et al. (2010) 



This paper introduces the first sample of optically- 
conf irmed SZE clusters from the ACT 2008 148 GHz data 
(see|High et al. 2010, for the recent SPT cluster follow-up 
study) and, for the first time, characterizes observation- 
ally the purity of an SZE survey. For this purpose we se- 
lected candidate galaxy clusters from ACT maps covering 
455 square degrees (a comoving volume of 2.46 Gpc^ at 
z < 1.2) using a rigorous and uniformly controlled proce- 
dure (see Marriage et al. 2010). Subsequently, during the 
2009B semester, we imaged candidates in the bands ^, r, 
and i on 4-m class telescopes to confirm the presence of 
a brightest cluster galaxy (BCG) and an accompanying 
red sequence of cluster members. 

In two companion papers we describe in detail the 
SZE cluster selection from the ACT data (Marriage 
et al. 2010) and cosmological constraints from the sam- 
ple presented here (Sehgal et al. 2010). Throughout 
this paper we assume a flat cosmology with Hq = 
lOOh km s~^ Mpc~^, h = 0.7, and matter density 
= 0.3. The reported magnitudes are always in 
the AB system. Cluster masses are quoted as M200, 
which corresponds to integrating the cluster mass within 
a spherical volume such that the density is 200 times the 
average density of the Universe at the cluster's redshift. 

2. OBSERVATIONAL SETUP 

Both ACT and SPT aim to provide a unique sample of 
massive clusters of galaxies, in principle selected by mass, 
nearly independently of redshift and over a large area 
of the sky covering a cosmologically significant volume. 
Utilizing these cluster samples to constrain cosmology 
and cluster physics requires empirical characterization of 
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the SZE selection function. 

As the first step in this process we measure the "pu- 
rity" of the ACT sample, i.e., the fraction of real clusters 
as a function of the SZE detection significance, using 
a signal-to-noise limited sample of cluster candidates ex- 
tracted from the 148 GHz maps. We search for an optical 
cluster associated with each candidate's SZE decrement. 
This is relatively straightforward, since for a wide range 
of mass limits in a ACDM cosmology 90% of clusters will 
lie below z ^ 0.8 and are therefore accessible in mod- 
est exposures on 4-m class telescopes. To illustrate this 
point, in Figure [l] we show the expected number of clus- 
ters of galaxies more massive than 4, 6 and 8 x IQ^ ^Mq as 
a function o f redshift from the mass function of ITinkerl 
et al. ( |2QQ8 ). This range in mass limit covers the cur- 
rent relevant range for ACT. A key requirement of our 
observing plan is efficient multi-band optical imaging of 
ACT cluster candidates. Unfortunately, the current gen- 
eration of large and wide format CCD mosaic cameras 
are not ideally suited for this task as their long read-out 
times (typically between 100 — 200 s) often become com- 
parable to the required exposure times to find L* (the 
characteristic Schechter luminosity) ellipticals ai z 1^ 
resulting in significant time overheads. For our observa- 
tions, we decided instead to trade the large field-of-view 
for smaller but faster readout CCD cameras (around 10 s) 
that allow us to move rapidly from target to target and 
secure multi-band imaging for the central 5' of the clus- 
ters, sufficient for the identification process because of 
the accurate astrometry of the mm-band data (see be- 
low). 

The number of nights required for this measurement is 
set by several factors, but the basic requirement is un- 
ambiguously detecting all optical clusters out to a red- 
shift of z ^ 0.8 (as shown in Figure |T]), where we can 
detect a high fraction (c^90%) of aumassive clusters 
even for a conservative mass threshold of 4 x 10^^ M©. 
The optical cluster detection requires seeing a cluster red 
sequence (i.e., early- type galaxies with luminosities less 
than L*). Our exposure times are therefore determined 
by the observed z— band magnitudes of low-luminosity 
(c^ 0.4L*) red-sequence galaxies at z 0.8, so that we 
can confidently identify clusters using the overdensity of 
red galaxies. In addition to being dominated by a red se- 
quence of galaxies, massive clusters exhibit another dis- 
tinctive signature, the presence of a BCG (typically with 
luminosity of a few L*). Our observation plan is hence 
designed to detect both of these features of real clusters 
(BCG and red-sequence) for 90% of the SZE candidates 
(i.e., z < 0.8). For most of the remaining clusters at 
z > 0.8 (estimated to be only 10% of the total number) 
we will be able to detect the presence of a BCG as well 
as some of the brightest cluster members. 

3. ANALYSIS 

3.1. Selection of ACT SZE Cluster Candidates 

ACT is a dedicated arcminute resolution CMB exper- 
iment consisting of a 6-m off-axis Gregorian telescope 
located at 5190m on Cerro Toco, near the Chajnantor 
Plateau, in the Atacama Desert in northern Chile (see 



from the ACT 2008 southern survey usi ng the Millime- 
ter Bolometer Array Camera (MB AC; [Niemack et al. 



2008|) on the ACT focal plane and with band centers 
Til 



arT48 GHz, 218 GHz and 277 GHz. The 2008 southern 
survey area is 9° wide, centered on declination —53.5°, 
and extends from right ascension 19^ to 24^ and 00^ 
to 07^30"^. A total of 850 hours of data with an aver- 
age of 680 working detectors at 148 GHz (yielding a data 
volume of 3200 GB) are used to make the microwave sky 
map. Relative detector pointings are established through 
observations of Saturn which has a high enough signal-to- 
noise ratio to allow a pointing fit for each detector. Ab- 
solute pointings for the survey are iteratively solved by 
comparing ACT radio source detections. The resulting 
rms uncert ainty in the reported positions is 5^' (Fowler 
[iFiI]|2QTQl). 

The data are calibrated with 6% precision to Uranus 
by extrapol ating the lower frequ ency WMAP7 Uranus 
calibration (Weiland et al. 2010). A highly-parallelized 
preconditioned conjugate gradient code is used to simul- 
taneously solve for the maximum likelihood millimeter- 
wave map and correlated noise (e.g., from large-scale at- 
mospheric emission). 

Cluster candidates were detected in a 148 GHz, 455 
square-degree submap from the ACT 2008 observing sea- 
son. The submap lies between right ascensions 00 12"^ 
and 07^08"" and declinations -56° 11' and -49°00'. The 
median noise level in this submap in terms of the decre- 
ment amplitude is 31/iK. After masking bright sources 
from the ACT 148 GHz source catalog, the map is match- 
filtered in the Fourier domain using the isothermal /3- 
model {(3 = 2/3) convolved with the ACT 148 GHz beam 
and models for the ACT noise, primary CMB fluctua- 
tions, and u ndetected sou rces. The form o f the filter 
dMelin et al.][2Q06, Haehnelt fc Tegmark|[T996, ) is 



Be^ {k) I Tother(A:) 



/ Be^ {k') I Tother(A:0 BeSk')dk' 



(1) 



Fowler et al.|2007[ [Hmcks et al.|[2QQ9l [Fowler et aIl[2Qro 



Swetz et al. |2010[ for a complete description). The SZE 
cluster selection tor this work is based on 148 GHz data 



where Bo^(k') is the Fourier transform of the beam- 
convolved /3-model with critical radius 0^^ and Tother(^0 
is the transform of the noise obtained from a combina- 
tion of ACT difference maps, the WMAP5 primary CMB 
signal, and the undetected sources (and c lusters) based 
on the amplitudes of Fowler et al. ( |2010 ). The map is 
filtered using /3-models with critical radii from 0.25' to 
4.0' in 0.25' steps. For a given filter, the SNR of a cluster 
detection is defined as the ratio of decrement amplitude 
to total rms. The reported significance of a cluster de- 
tection corresponds to the filter which maximizes this 
SNR. Furthermore, any detection that was not detected 
at greater than SNR>3.5 in at least four conjoined pixels 
was discarded. 

During 2009B the ACT data reduction pipeline and 
cluster extraction algorithms evolved such that every- 
thing from pointing to data selection to the matched 
filter for the final candidate list changed from one ob- 
serving run to the next. Thus we observed more SZE 
candidates than used in the purity study (although we 
report all confirmed cluster detections regardless of its 
origin). Despite this inefficiency the optical study com- 
pletely vetted a SNR-ordered list of candidates selected 
in the manner described above. A final caveat is that af- 
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Table 1 

2009B Approved Observing Runs 



Run 


Tel. 


Camera 


Dates 


Nights 


TAG 


1 


SOAR 


SOI 


2009, Aug 22-23 


0/2 


GNTAG 


2 


NTT 


EFOSC 


2009, Oct 22-26 


3/4 


ESO/Ghile 


3 


SOAR 


SOI 


2009, Dec 9-12 


4/4 


NOAO 



Note. — For run 1 we applied for SOAR time but were reas- 
signed by the GNTAG to the GTIO/Blanco 4-m telescope due to 
re-aluminization of the SOAR mirror. 

ter the optical observing season ended, the ACT pipeline 
continued to evolve such that the finalized candidate list 
tested here is not identical to the list presented in Mar- 
riage et al. (2010). While the SNR of the high signifi- 
cance candidates remains similar between this study and 
Marriage et al. (2010), the fainter candidates differ in 
SNR with some of the faintest candidates falling below 
the SNR threshold of the Marriage et al. list. This is ex- 
pected: different map-making configurations weight the 
data in different ways, and a faint candidate which is 
decrement-boosted by noise in one realization of map 
may not be boosted in the next. 

3.2. 2009B Optical Observing Campaign 

During the 2009B semester we were awarded a total of 
10 dark nights for this project in three observing runs on 
4-m class telescopes obtained through different Telescope 
Allocation Committees (TACs). Unfortunately our first 
run (2 nights), which occurred during the Chilean winter, 
was completely lost due to snow and no data were ob- 
tained (09B-1007, PLInfante). The next two observing 
runs had optimal conditions and enabled the successful 
completion of our program. In Table [l] we summarize the 
telescopes, instruments, dates, number of good (photo- 
metric) nights, and the TACs that awarded the time. 

In the run during October 22-26, 2009 we used the ESQ 
Faint Object Spectrograph and Camera (EFOSC) on the 
3.6-m NTT on La Silla (084.A-0619, PLBarrientos) and 
observed SZE candidates using the Gunn gri filter, with 
exposure times of 270 s (3 x 90 s), 360 s (3 x 120 s) and 
1100 s (4 X 275 s) in g, r and i respectively. EFOSC 
consists of a 2048 x 2048 pixel thinned, AR coated 
CCD#40 detector projecting a 4.1' x 4.1' field-of-view 
at 0.2408'7pixel (2x2 binning) on the sky. For the 
three first nights of this run atmospheric conditions were 
photometric with seeing 0.''7 — 1" . The last and fourth 
night of the run was lost due to clouds and no data were 
taken. Throughout the run we also observed photomet- 
ric standards from the Southern Hemisphere Standards 
Stars Catalog (Smith et al. 2007) to calibrate our data. 

For our final run, we used the SOAR Optical Imager 
(SOI) on the 4.1-m SOAR Telescope on Cerro Pachon 
(09B-0355, PI: Menanteau) for four nights from Decem- 
ber 9-12, 2009. SOI is a bent-Cassegrain mounted optical 
imager using a mini-mosaic of two E2V 2048 x 4096 pixel 
CCDs to cover a 5.26' x 5.26' field-of-view at a scale of 
0."154/pixel (2x2 binning). For our cluster search we 
used the same exposure times as in the NTT using the 
SDSS gri filter set. For three cluster candidates pre- 
sumed to be at high redshift we obtained observations 
twice as long as the nominal exposure times and addi- 
tional deep 2200s (8 x 275s) z-band observations. The 
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SOAR/SOI 
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Figure 2. The recovery fraction of simulated elliptical galaxies, 
using a de Vaucouleurs profile, as a function of the observed i- 
band magnitude for the NTT (dashed line) and SOAR (solid line) 
observations. We show for reference the 50% and 80% recovery 
fraction levels. 

run conditions were optimal and ah four nights were pho- 
tometric, with seeing in the range 0/^5 — \" . We ob- 
served photometric standard stars using the same set and 
method as we did for the earher NTT run. 

3.3. Optical Data Reduction 

For the optical analysis we take advantage of an ex- 
isting imaging data pipeline that takes us from raw 
data to galaxy catalogs. We use a modified version 
of the Python- based Rutgers Southern Cos mology im- 
age pipeline (Menanteau et al. 2009 2010) to process 
the data from the NTT and SOAR runs. The ini- 
tial standard CCD processing steps are performed using 
IRAFpVccdproc tasks via the STScI/PyRAfj^ inter- 
face. The fast read-out time of the CCDs allowed us to 
accommodate multiple dithered exposures for each tar- 
geted cluster candidate. In our cluster observations we 
followed a 10^' off-center 3-, 3-, and 4-point dither pat- 
tern in the ^, r, and i bands, respectively, which helped 
to avoid saturated stars, remove chip defects and arti- 
facts, and cover the 10^' gap between CCDs in the SOI 
camera. Individual science frames were astrometrically 
recalibrated before stacking by matching sources with 
stars from the US Naval Observatory Catalog. Observed 



fields with p hotometric standard stars from the [Smith 
et al. (2007) catalog were processed in the same fashion 
and we obtained a photometric zero-point solution in AB 
magnitudes for each observing night. For each targeted 
cluster we associated all science frames in all filters and 
aligned and median combined them int o a co mmon pixel 
reference frame using SWarp (Bertin 2006| ) to a plate 
scale of 0.''2408/pixel for the NTT and 0.''1533/pixel 
for SOAR. Source detection and photometry f or the sci- 
ence catalogs were performed using SExtractor ( Bertin & 



, IRAF is distributed by the National Optical Astronomy Ob- 
servatory, w hich is operated by the AUHA under cooperative a^ree- 
ment with NiSi:^' 



~^^~PyRAi:^' is a product ot the iSpace I'elescope iScience Institute. 
^rhich IS operated by AUHA for 
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Figure 3. The observed i-band magnitudes of L*, 0.4L* and 4L* 
(BCG) early- type galaxies as a function of redshift. We use L* as 
defined Blanton et al. (2003) for the population of red galaxies at 
z = 0.1 and allow it to passively evolve with redshift. We show 
in different gray levels the 80% completeness limits for the NTT 
and SOAR observations. In other words, at magnitudes below the 
gray band we detect elliptical galaxies with at least 80% complete- 
ness. We also show the observed magnitudes for the three brightest 
galaxies in three high redshift clusters. 

|Arnouts'1996) in dual-image mode in which sources were 
identified on the z— band images using a 1.5<j detection 
threshold, while magnitudes were extracted at matching 
locations from all other bands. The computed magni- 
tudes in our catalogs were corrected for Galactic dust 
absorption in each observed band utilizing the i nfrared 
maps and C routines provided by Schlegel et al.| ([l998 ) . 
This step is particularly important tor obtaining unbi- 
ased colors, required for accurate photometric redshift s, 
as significant variations occur over the large area covered 
by the ACT 2008 survey region. 

The final step in our pipeline involves computing the 
photometric redshift and redshift probability distribu- 
tions, Pbpz{z)^ for each object using the spectral energy 
distribution (SED) based BPZ code (Benitez 2000) with 
a fiat prior. We compute photometric redshift s using the 
dust-corrected ^, r, i isophotal magnitudes, as defined by 
the 2-band detection and the BPZ se t of template spec- 
tra described in Benit ez et al. | ([2QQ4 | which is based on 
the templates from Coleman et al. l [T980p and [Kinney et 
g1 (Il996). This set consists orEl,^Ec, Scd, im, ^B3, 
and SB2 and represent the typical SEDs of elliptical, 
early /intermediate type spiral, late- type spiral, irregu- 
lar, and two types of starburst galaxies respectively. 

3.4. Galaxy Completeness and Recovery Fraction 

The exposure times for our observations were designed 
to confidently detect sub-L* early-type galaxies up to 
z ~ 0.8 and BCGs at z > 1.0 in order to confirm the 
presence of a galaxy cluster. In order to assess the depth 
reached by the observations as well as the differences be- 
tween telescopes, we performed Monte Carlo simulations 
to estimate the recovery fraction of early-type galaxies 
as a function of magnitude in each of our observing runs. 
We used the IRAF/artdata package to create simulated 
data sets for the NTT and SOAR observations, which 



were subsequently processed through our pipeline using 
the same detection and extraction procedures as for the 
science images. Within artdata we used the mkobjects 
ta sk to draw artifi cial elliptical galaxies on science fields 
as|de Vaucouleurs| profiles with the same noise properties 
Ve scaled the surface brightness of the 



galaxy profiles according to their magnitude for a range 
of assumed sizes (1 — 7 kpc, but converging to the typical 
image seeing at the faint limits). Using this procedure 
we generated 20 and 25 artificial elliptical galaxies for 
the NTT and SOAR datasets respectively over a wide 
range of magnitudes {19 < i < 27) sampled at small 
(Amag = 0-1) steps. We repeated this 10 times for each 
dataset at each step. At each magnitude step we matched 
the recovered galaxies in the simulations with the input 
positions from the artificial catalogs. In Figure[2]we show 
the recovered fraction of galaxies as a function of the in- 
put i-band magnitude for the NTT and SOAR dataset, 
as well as the 80% and 50% levels for reference. From this 
exercise we conclude that we reach 80% completeness at 
i = 23.5 and i = 24.3 for elliptical galaxies for the NTT 
and SOAR datasets respectively. We emphasize that we 
followed the same observing strategy (i.e. same filter set, 
dithering pattern and exposure times) for both runs and 
therefore the relatively large difference (0.8 mag) in the 
completeness limit is likely due to a combination of tele- 
scope size (3.6-m vs 4.1-m), CCD sensitivity and tele- 
scope throughput. The re-aluminization of the SOAR 
primary mirror in November 2009, just before our run, 
might have played a role in the higher sensitivity of the 
SOAR imaging. 

We can use the completeness limits estimated from our 
simulations to determine how far in redshift we can "see" 
massive clusters. For this, we compare the completeness 
limits of our observations to the expected and observed 
(i.e. known) apparent magnitudes of galaxies in clusters 
as a function of redshift. We estimated the expected ap- 
parent galaxy z-band magnitude as a function of redshift 
using L* as de fined for the population of red galaxies by 
Blanton et al.||2003) at z = 0.1 and allowing passive evo- 
lution according to a solar metallicity Bruzual & Cha riot] 
(2003) r = 1.0 Gyr burst model formed at Zf = 5. We 



show this in Figure |3]for a range of luminosities (0.41/ 
L* and 41/*) aimed at representing the cluster members 
from the faint ones to the BCG. We also show as different 
gray levels the 80% completeness as determined by the 
simulations for the NTT and SOAR run, which provides 
evidence that we can reach typical L* galaxies to z ~ 0.8 
for both datasets (and z ~ 1 for the SOAR data). We 
also compare the depth of our imaging campaign to the 
observed i-band apparent magnitudes of galaxies in three 
massive high redshift {z > 1.2) clusters for which we 
could find suitable data in the literature for comparison. 
In Figure [3] we plot the observed magnitude of the three 
brightest galaxies for the clusters RDCS 1252.9-2927 



(iBlakeslee et al. 



2003 



Rosati et al.ll2004t at z = 1.237 

05 



XMMU J2235.T 

at z = 1.393 and XMMXCS J 2215.^ 



Rosati et al. 

. ^ ^ -1738 ( IStan- 

Iford et al.||2006l [Hilto n et al.|[2QQ9| at z = 1.457"^ 

conclude that, with the possible exception of a cluster 
like XMMXCS J2215.9-1738 observed during the NTT 
run, we comfortably detect luminous galaxies in clusters 
at high redshift. This provides further evidence that our 
observation's integration times allow us to confidently 
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Table 2 

ACT 2008 clusters 



ACT Descriptor 


Vf \ ( To^^^^ 
il.A. (^JzUUUj 


T^r^/^ ( To^^^^ 
uec. (^JZUUUJ 


xieasniit 


olN IX 


Alt iName 


AC i-CL 


J0145-5301 


01:45:03.6 


-53:01:23.4 


0.118 ^ 


4.7 (4.0) 


Abell 2941 


ACT-CL 


J0641-4949 


06:41:37.8 


-49:46:55.0 


0.146 ^ 


4.9 (4.9) 


Abell 3402 


ACT-CL 


J0645-5413 


06:45:29.5 


-54:13:37.0 


0.167 ^ 


7.1 (7.1) 


Abell 3404 


ACT-CL 


J0638-5358 


06:38:49.4 


-53:58:40.8 


0.222 ^ 


10.6 (10.0) 


Abell S0592 


ACT-CL 


J0516-5430 


05:16:37.4 


-54:30:01.5 


0.294 ^ 


5.2 (4.7) 


Abell S0520/SPT-CL J0516- 


ACT-CL 


J0658-5557 


06:58:33.1 


-55:57:07.2 


0.296 ^ 


11.6 (11.5) 


1E0657-56 (Bullet) 


ACT-CL 


J0245-5302 


02:45:35.8 


-53:02:16.8 


0.300 


8.3 (9.1) 


Abell S0295 


ACT-CL 


J0217-5245 


02:17:12.6 


-52:44:49.0 


0.343 ^ 


4.5 (4.1) 


RXC J0217. 2-5244 


ACT-CL 


J0237-4939 


02-37-01 7 


-4Q-38-10 




4.9 (3.9) 




ACT-CL 


J0707-5522 


07:07:04.7 


-55:23:08.5 


0.43 ± 0.06 


4 2 f ) 




ACT-CL 




02:35:45.3 


-51:21:05.2 


0.43 ±0.07 


5.7 (6.2) 




ACT-CL 


J0330-5227 


03:30:56.8 


-52:28:13.7 


0.440 g 


7.4 (6.1) 


Abell 3128(NE) 


ACT-CL 


J0509-5341 


05:09:21.4 


-53:42:12.3 


0.461 h 


4.4 (4.8) 


SPT-CL J0509-5342 


ACT-CL 


J0304-4921 


03:04:16.0 


-49:21:26.3 


0.47 ±0.05 


5.0 (3.9) 




ACT-CL 


J0215-5212 


02:15:12.3 


-52:12:25.3 


0.51 ±0.05 


4.8 (4.9) 




ACT-CL 


J0438-5419 


04:38:17.7 


-54:19:20.7 


0.54 ±0.05 


8.8 (8.0) 




ACT-CL 


J0346-5438 


03:46:55.5 


-54:38:54.8 


0.55 ±0.05 


4.4 (4.4) 




ACT-CL 


J0232-5257 


02:32:46.2 


-52:57:50.0 


0.59 ±0.07 


5.2 (4.7) 




ACT-CL 


J0559-5249 


05:59:43.2 


-52:49:27.1 


0.611 ' 


5.1 (5.1) 


SPT-CL J0559-5249 


ACT-CL 


J0616-5227 


06:16:34.2 


-52:27:13.3 


0.71 ±0.10 


6.3 (5.9) 




ACT-CL 


J0102-4915 


01:02:52.5 


-49:14:58.0 


0.75 ±0.04 


8.8 (9.0) 




ACT-CL 


J0528-5259 


05:28:05.3 


-52:59:52.8 


0.768 h 


4.7 (...) 


SPT-CL J0528-5300 


ACT-CL 


J0546-5345 


05:46:37.7 


-53:45:31.1 


1.066 h 


7.2 (6.5) 


SPT-CL J0546-5345 



Note. — R.A. and Dec. positions denote the BCG position in the optical images of the cluster. The SZE position was used to construct 
the ACT descriptor identifiers. Clusters with uncertainty estimates on their redshifts are those systems for which only photometric redshifts 
are available from the NTT and SOAR gri imaging. The values of SNR in parentheses are the current values from Marriage et al. (2010), 
while the others are the ones used in this study 
^ spec-z from de Grandi et al.|_|1999|) 
^ spec-z from Jones et al. (200^ 
^ spec-z from Guzzo et al.„(,19g 
^ spec-z from Tucker et al.J (|1998| > 
^ spec-z from 'Edge et al. 7-^994') 
^ spec-z from Bohringer et al. (2004]) 
s spec-z from Werner et al. (20U7j) 
^ spec-z from Infante et al. ( 201o| ), [Brodwin et al. 

" " ~'m 



(20101 



^ spec-z from High e t al.| ( |2010l > 

identify real SZE clusters up to z ^ 1.2 — 1.4, below this 
redshift according to ACDM model > 99% of all c lust ers 
more massive 8 x 1{)^^M(?) should reside. In Section 4.7 we 



further discuss the possibility of missing z > 1.2 massive 
clusters in our observations. 

3.5. Additional Imaging Data 

A number of the SZE cluster candidates were associ- 
ated with previously known clusters at lower redshifts. 
For a few of these we obtained new NTT/SOAR gri ob- 
servations in order to calibrate our photometric redshift 
and to test our identification procedure; however, for the 
most part these were avoided and we relied on public 
archival imaging and on data from other existing pro- 
grams to compare with the SZE sources. 

AS0295 and AS0592 are incl uded amon^ the low red- 
shift SZE clusters reported by [Hincks et al. (2009j) and 
are the focus of a separate observmg program by our 
group (09B-0389, PLHughes) aimed at obtaining weak 
lensing masses using deep V and R observations. These 
were obtained using the MOSAIC camera on the Blanco 
4-m telescope on January 9, 2010 under photometric con- 
ditions. 

1E0657— 56 (the Bullet Cluster) has been extensively 
observed at different wavelengths. Here we show the 
optical observations of the cluster central region taken 



with the Advanced Camera for Surveys aboard the Hub- 
ble Space Telescope (HST) (GO: 10200, PI: Jones and 
GO:10863, PLGonzalez) using the F606W (V), F814W 
(/) and F850LP {z) filters which were taken from the 
HST archiv4!!l 

For RXCJ0217.2-5244 we used our own pipeline to 
process raw V (720s) and R (1200s) images obtained from 
the ESO archiv^ (70.A-0074-A, PLEdge) utilizing the 
SUSI camera on the NTT on Oct 7 and 12, 2002. 

A fraction of the area covered by ACT overlaps with 
the 41 deg^ griz imaging from the CTIO 4-m telescope 
MOSAIC camera on t he 05hr field from th e Southern 
Cosmology Survey (Menanteau et al. 2010). We used 



these data to search for optical identitication of SZE clus- 
ter candidates in the region to avoid re-targeting with the 
NTT and SOAR telescopes. 



4. RESULTS 



4.1. 



The ACT SZE Cluster Sample 

Our analysis has resulted in a sample of 23 optically- 
confirmed SZE clusters selected from 455 square-degrees 
of ACT data between 0.1 < z < 1.1. Nine of these 
clusters, such as 1E0657— 56 (the Bullet Cluster), Abell 

http://archive.stsci.edu/ 
http://archive.eso.org/ 
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Figure 4. Composite color images for ACT-CL J0145-5301 (Abell 2941), ACT-CL J0641-4949 (Abell 3402), ACT-CL J0645-5413 (Abell 
3404) and ACT-CL J0638-5358 (Abell S00592). The horizontal bar shows the scale of the images, where north is up and east is left. White 
contours show the 148 GHz SZE maps with the minimum and maximum levels, in /iK, displayed between brackets. 

S0592, and Abell 3404 are well known and appear at low 
redshift [z < 0.4); and four systems in our sample have 
been reported as SZE-selected clusters by SPT in their 
area overlapping with the ACT coverage. Three of these 
systems (ACT-CL J0509-5341, ACT-CL J0528-5259, 
and ACT-CL J0546— 53 45) were orifflnally repor ted as 
SZE-selected clusters by Staniszewski et al. (2009), with 
optical and X-ray properties subsequently studied by us 
dMenanteau fc Hughes|20Q 9 [), and AC T J0559-5249 has 
been recently reported by High et al. (2010). About two- 
thirds of the clusters we present here (14), however, are 
new SZE-discovered systems (4 of which overlap with 



the SPT cluster sample of High et al. 2010) previously 
unrecognized in the optical or X-ray bands; ten of these 
are entirely new discoveries by ACT, with photometric 
redshifts ranging from z = 0.4 to z = 0.8. All of the new 
clusters are optically quit e rich and a number show strong 
lensing arcs (see section 4.2). We present the full list 
of clusters, positions and redshifts (photometric in some 
cases) in Table |2j In Figures 4pQ we show the composite 
color optical images of all clusters in our sample with 
the ACT 148 GHz contours (six equally-spaced between 
the minimum and maximum levels in /iK shown between 
brackets) overlaid on the images. 
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ACT-CL J0245-5302, z=0.300 
(AS0295) 



ACT-CL J0217-5245, z^O.343 
(RXCJ0217. 2-5244) 




Figure 5. Composite color images for ACT-CL J0516-5430 (Abell S0520), ACT-CL J0658-5557 (Bullet), ACT-CL J0245-5302 (Abell 
S0295) and ACT-CL J0237-4939 (RXC J0217.2-5244). The horizontal bar shows the scale of the images, where north is up and east is left. 
White contours show the 148 GHz SZE maps with the minimum and maximum levels, in /iK, displayed between brackets. 



In the following sections we provide detailed informa- 
tion on individual clusters in the sample in increasing 
redshift order. 

4.1.1. ACT-CL J0145-5301 (Abell 2941) 

Abell 2941 is the lowest redshift cluster in o ur sample 
and it was originally reported by Sersic ( 1974 ) as a com- 
pact cluster in his list of peculiar galaxies, groups and 
clusters based on the Maksutov plate collection of the 
Universidad de C hile's Observatorio Astronomico. It was 
later included by "Abell et al. ( 1989 ) as an optically rich 
cluster, with R=2 according to the Abell ( 1958} classifi- 



cation. Abell 2941 has also be en identified a s a ROSAT 
Ah Sky Survey (RASS) source ( |Voges et al.fl 999) and X- 
ray cluster RXCJQ145-53QQ dz = 0.1183 by |de Grandi] 
etaD(IT999l. 



4.1.2. ACT-CL JO64I-4949 (Abell 3402) 

Abell 34 02 is included in t he catalog of rich optical 
clusters of [Abell et al. (1989) with R=0 richness. We 
identify its BCG (by exammmg the DSS imag es) as the 
2MASS sou rce 2MASS J06413783-4946548 ( |Skrutskie| 
et al.|[2QQ6| ). The cluster has not been targeted spectro- 
scopically, however, the redshift of the BCG {z = 0.146) 
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Figure 6. Composite color images for ACT-CL J0237-4939, ACT-CL J0707-5522, ACT-CL J0235-5121 and ACT-CL J0330-5227 (Abell 
3128 North-East). The horizontal bar shows the scale of the images, where north is up and east is left. White contours show the 148 GHz 
SZE maps with the minimum and maximum levels, in /xK, displayed between brackets. 



ha s been reported as part of the 6dF survey (Jones et 
aLl|2009|). It is worth mentioning that despite its low 



redshitt, Abell 3402 has not been reported as an X-ray 
source by any survey yet (including the RASS). More- 
over, the relatively large distance between the SZE cen- 
troid and the cluster BCG (1.8^) may suggest that this 
is potentially a spurious association. Note that the BCG 
falls right at the edge of the field of view shown in Fig. 4 
(upper right panel). Further spectroscopy and X-ray ob- 
servations of the cluster will provide definite information 
to determine whether this is a bona fide association with 
the SZE detection or a case of a "cluster behind a cluster" 



like ACT-CL J0330-5227 (see section 4.1.8) 



4.1.3. ACT-CL JO645-5413 (Abell 3404) 

Abell 3404 was first reported by |Abe ll et ah] ( |1989[ ) as 
an optical cluster with R=l richness cla ss ancTit was als o 
detected as a RASS bright source (Voges et al. 19991). 
The cluster redshi ft, z = 0.1661, was reported by |ae| 
Grandi et al. (1999) and the system is also know b y its 
R El^'LEX designation RXCJ0645.4-5413 (Bohringeret] 
aI][2QQ4l ). 



4.1.4. ACT-CL J0638-5358 (Abell S0592) 
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Figure 7. Composite color images for ACT-CL J0509-5341 (SPT-CL J0509-5342), ACT-CL J0304-4921 , ACT-CL J0215-5212 and ACT- 
CL J0438-5419. The horizontal bar shows the scale of the images, where north is up and east is left. White contours show the 148 GHz 
SZE maps with the minimum and maximum levels, in fiK, displayed between brackets. 



Abell S0592 was first reported by l Abell et al.| ( |1989D 
in tiieir supplementary table of southern clusters. THese 
were clusters not rich enough or too distant to satisfy the 
criteria to be included in the rich nearby main cluster cat- 
alog. It was also later reported as a ROS AT hnght source 
in the RASS ( Voges et al. 1999) and as the RE FLEX X- 
ray cluster, R XC J 0638.7^5 8 (Bohr inger et al. 2004) 
at z = 0.2216 dde Grandi eFal.,,1999^. Abeh S0592 was 
first rep orted as an ^ZE sour ce by ACl' (ACT-CL J0638- 



Abell et al. 



Abell S05 2 is a n optically-rich cluster reported by 
(1989) in their supplementary c atalog and 

the 



it is located at 2; = 0.294 (IGuzzo et al.|[T999i). It has 



been reported as a ROS AT bright source in the RASS 



fVoges et al.l 11999) and it is kno wn by its REFLEX 



5358) in Hincks et al. (2009). 



designation RXCJ0516.6-5430 (Bohringer et al ||2004[) 
This cluster was also detec ted as an SZE cluster by SFT 
( Staniszewski et al. 2009) and the results from several 
mass estimates (optical, X-ray and weak lensing) sug- 
gest a value of M > 5 x 10^"^ M^:, (Zhang et al.||2006' 



4.1.5. ACT-CL J0516-5430 (Abell S0520) 



Menanteau fc Hughes|[2QQ9l [Mclnnes et al.„2009| ): 
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Figure 8. Composite color images for ACT-CL J0346-5438, ACT-CL J0232-5257, ACT-CL J0559-5249 ( SPT-CL J0559-5249) and ACT- 
CL J0616-5227. The horizontal bar shows the scale of the images, where north is up and east is left. White contours show the 148 GHz 
SZE maps with the minimum and maximum levels, in /xK, displayed between brackets. 

4.1.6. ACT-CL J0245-5302 (Ahell S0295) 



et al. 



re ported as an ACT SZ E source (ACT-CL J0638-5358) 
by Hincks et al. ( |2QQ9|) . ASCA observations yielded an 
average tem perature oi /cT = 6.72 ± 1.09 keV ( Fukazawal 
1994|), was originally selected optically (Abell e| [ittL]^04|. Figure |5] shows our deep, multicolor optical 



The galaxy cluster Abeh S0295 at z = 0.300 OEdge 



al.l|l989') and was also detected in the RASS (Voges ef 



al.|l U99 ) . There has been a reported discovery of a giant 
strong- lensing arc near the brightest cluster galaxy (Edge 
|et al. 199 4), which we confirm in our imaging (see lower 
left panel of Figure |5|. There were also previous (unsuc- 
cessful) attempts at de tecting the SZE effect at 1.2-mm 
and 2-mm using SEST ( [Andreani et al.|1996D . The ACT 
detection of this cluster is highly signiticant (it is ranked 
in the top 5 of the sample of 23). Abell S0295 was first 



imaging of this cluster from the CTIO 4-m telescope in 
January 2010. The distribution of galaxies is highly elon- 
gated and matches well the shape of the SZE contours; 
both suggest the idea that this is a merging system. 

4.1.7. ACT-CL J0217-5245 (RXC J0217.2-52U) 

RXC J0217.2-5244, at z = 0.343, was first identified 
as an X-ray cluster by the REFLEX survey (Bohringer et| 
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Figure 9. Composite color images for ACT-CL J0102-4915 and ACT-CL J0528-5259 (SPT-CL J0528-5300) The horizontal bar shows 
Fhe scale of the images, where north is u p and east is lett. Wiiite contours sJiow tJie 14^ UhLz bZiiil maps witJi tiie minimum and maximum 
levels, in /iK, displayed between brackets. 




Figure 10. Composite color images for ACT-CL J0546-5345 
(SPT-CL J0546-5345). The horizontal bar shows the scale of the 
images, where north is up and east is left. White contours show 
the 148 GHz SZE maps with the minimum and maximum levels, 
in fiK, displayed between brackets. 



Abell 3128 North-East (NE) is a great example of the 
power of the SZE to detect clusters regar dless of red- 
shift a nd only limited by mass. Until , Werner et al. 
(2007), using XMM-Newton data, revealed a more dis- 
tant [z = 0.440) and more massive cluster superposed 
on the northeastern component of Abell 3128, it was be- 
lieved to be part of the Horologium- Reticulum superclus- 
ter at z = 0.06, although its X-ray morphology is clearly 
double peaked wit h the two compo n ents separated on the 
sky by some 12^ (|Rose et"aL]|2002D . [Hincks et aLj ( |2009 ) 
confirmed thel WerneT 3t al .|(|2007p results that AbeOT^ 
(NE) is indeed a separate higher redshift cluster, and a 
case of finding a cluster behind a cluster. We do not 
detect with ACT a significant decrement from the lower 
redshift southwestern component. We show our optical 
EFOSC/NTT imaging of the cluster in Figure (6] There 
is a strong lensing arc a few arcsecond southwest of the 
central brightest cluster galaxy. 

4.1.9. ACT-CL J0658-5557 (1ES0657-558) 

1ES0657— 558 is the famous "Bullet" cluster, which we 
detect with the highest significance of any cluster in the 
ACT sample. Previous detections of the mm-band SZE 
signal from this clus ter have bee n reported by ACBAR 



dComez et al."2004 l), APEX- SZ ( |Halverson et al. 
and SP'i' (^lagge et al.||2010D . 



2009) 



al 



2004). Close comparison of the SUSI/NTT archival 



imagmg with the ACT 148GHz contours (see Figure [5| 
shows a double-peaked morphology in the SZ maps that 
matches the position of the cluster BCG and the second 
brightest elliptical galaxy. 

4.1.8. ACT-CL J0330-5227 (Ahell 3128 North-East) 



4.1.10. Overlap with the SPT sample 

Four clusters in our sample overlap with SZE discov- 
ered systems previously reported by SPT. This presents 
an opportunity to compare observed clusters properties 
using different observing strategies and techniques (al- 
beit sometimes the same archival data). 

ACT-CL J0509-5345, ACT-CL J054 7-5345 and 
ACT- CL J0528-5399 were first reported by 'Staniszewskij 
et al. (2009j) and their physical properties (photometric 
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Figure 11. Close-up monochrome images of examples of strong lensing arcs in three clusters in the ACT sample. 



redshifts, luminosities and m ass estimates) reported in 
Menanteau fc Huglie's| ([2009'). They have been the fo- 
cus of intense multi-wavelength observations to estab- 
lish their redshifts and estimate their masses using dif- 
ferent techniques. Our team has secured multi-object 
spectroscopic redshifts for these systems using F0RS2 
on the VLT during 2009B (084.A-0577, PLInfante) as 
part of a dedicated program aimed to calibrate the 
observable-cluster mass scaling relations of SZE clus- 
ters at high redshift (see Infante et al. 2010). Two of 
these systems, J0509-5345 {z = 0.461) and J0528-5399 
{z = 0.768), have spectros copic redshifts con sistent with 
those recently reported by [High et al. (2010). For ACT- 
CL J0547— 5345 we have secured the highest spectro- 
scopic redshift, z = 1.0678 for an SZE-sel ected clusters 
until now (Infante et al. 2010, Brodwin et al. 2010). 

The spectroscopic redshift for ACT-CL J0559-5249, 
(SPT-CL J0559- 5249), z = 0.611 has been reported by 



High et al.|2010| and provides a consistency check for our 
photometric redshift for this system, z = 0.56 ± 0.06. 

4.2. Gravitational lensing features 

A number of the clusters in our sample show gravi- 
tational lensing features, indicative of massive systems. 
The excellent seeing conditions of the 2009B runs used 
for this study, enable us to unambiguously determine 
that 6 of the 23 clusters show strong lensing arcs. 
Three of these have already been reported in the liter- 
ature: ACT-CL J0658-5557 (lE0657-56/Bullet), ACT- 
CL J0245-5302 (AS0295), and ACT-CL J0546-5345. 
Arcs in the other three are reported here for the first 
time: ACT-CL J0330-5227 (Abeh 3128-NE), ACT- 
CL J0509-5341, and ACT-CL J0304-4921. In Figure [11 
we show close-up monochromatic images of the clusters 
with newly identified strong lensing arcs in the sample. 

4.3. Photometric Redshifts of the Clusters 

Ten systems in our sample of ACT SZE clusters are 
completely new systems, for which the only optical data 
available are from our NTT and SOAR imaging collected 
during 2009B. The main objective of our observations 
was to confirm SZE clusters via optical identification as 
real systems; and in the case that observations were car- 
ried under photometric conditions — which was true for 
both NTT and SOAR runs — to compute crude photo- 
metric redshift estimates using the gri (and z in some 



cases) imaging for the observed clusters. Therefore we 
can provide only photometric redshift estimates for these 
new systems (see Table [2|. A small subsample of the 
targeted clusters (9) has spectroscopic information avail- 
able, which we can use to perform some basic compari- 
son between the photometric and spectroscopic redshifts. 
Unfortunately, the imaging of this subsample is quite het- 
erogeneous and comes from different datasets (4 from 
CTIO 4-m, 3 from NTT and 2 from SOAR) making it 
difficult to characterize the photometric redshift errors 
for the sample. In Figure [12] we show the comparison 
between spectroscopic and photometric redshift for the 
mean cluster value (top panel) and the BCG (bottom 
panel) for the subsample of nine clusters. We calculate 
the mean photometric redshift for each cluster, ^^duster, 
wit h a pr ocedure similat to that used in Menanteau et| 
al. ( |201 0|). This consists of iteratively selecting cluster 
member galaxies using a 3a median sigma-clipping algo- 
rithm within a local color-magnitude relation (using all 
available colors) defined by galaxies that: a) were pho- 
tometrically classified as E or E/SOs according to their 
BPZ-fitted SED type, b) are within a projected radius 
of 500 kpc from the BCG and c) are within the redshift 
interval \z — ^^ciusterl = = 0.05. The uncertainty on 
this redshift estimate was taken to be the weighted rms 
of the individual galaxies chosen as members. The uncer- 
tainties associated with the BCG photometric redshifts 
were estimated from the BCG Pbpz(^) function. These 
probability functions are non- Gaussian as revealed by the 
asymmetric error bars in Figure 12 (lower panel). 

In Table [3] we provide the mean cluster and BCG pho- 
tometric redshift estimates and associated 68% C.L. un- 
certanties. Despite the relatively large errors on the pho- 
tometric redshifts (driven by the limited number of bands 
used) they are still useful for providing cosmological con- 
straints (see Sehgal et al. 2010). 

4.4. False Positive Rate 

We expect fewer than 0.010 false positives based on tar- 
geting 49 regions of the sky and assuming an expected 
cluster number density of 0.06 deg~^ and a positional 
offset between the BCG and the SZ target centroid of 2^ 
This assumes a random association with a cluster at or 
above the mass threshold of our sample. We would ex- 
pect to find one random association for a cluster density 
of 5.96 deg~^, which corresponds to a mass threshold of 
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Table 3 

Cluster Photometric Redshifts 



ACT Descriptor 


^cluster 


ZBCG 


Zspec 


Tel. 


Filters 


ACT-CL J0232-5257 


0.59 ± 0.07 


'-'•^^-0.63 




EFOSC/NTT 


9,r,i 


ACT-CL J0304-4921 


0.47 ±0.05 


^ , ^J-D 41 
^•^^-0.51 




EFOSC/NTT 


9,r,i 


ACT-CL J0330-5227 


0.52 ± 0.07 




0.440 


EFOSC/NTT 




ACT-CL J0438-5419 


0.54 ±0.05 


_ _ . _i_n 44 

"^•^^-0.56 




EFOSC/NTT 


9,r,i 


ACT-CL J0616-5227 


0.71 ±0.10 


^- ' -"--0.78 




EFOSC/NTT 


9,r,i 


ACT-CL J0641-4949 


0.18 ± 0.07 


^ _ ^-LO 1 4 


0.146 


EFOSC/NTT 


9,r,i 


ACT-CL J0645-5413 


0.19 ± 0.04 


^ . „-i-n 1 f! 

'-'•-^^-0.22 


0.167 


EFOSC/NTT 


9,r,i 


ACT-CL J0102-4915 


0.75 ± 0.04 


0.78l°-i« 




SOI/SOAR 


9,r,i,z 


ACT-CL J0145-5301 


0.12 ± 0.06 


^ ^ ^^-l-D 1 1 


0.118 


SOI/SOAR 


9,r,i 


ACT-CL J0215-5212 


0.51 ± 0.05 


n r-o+O.SO 

0-52Io.55 




SOI/SOAR 


9,r,i 


ACT-CL J0235-5121 


0.43 ± 0.06 


42+0-31 




SOI/SOAR 


9,r,i 


ACT-CL J0237-4939 


0.40 ± 0.05 


U.C)D_Q 43 




SOI/SOAR 


9,r,i 


ACT-CL J0346-5438 


0.55 ±0.05 


'-'•^^-0.64 




SOI/SOAR 


9,r,i 


ACT-CL J0559-5249 


0.56 ±0.06 


56+0-55 

'-'•^"-0.64 


0.611 


SOI/SOAR 


9,r,i,z 


ACT-CL J0707-5522 


0.43 ±0.07 


44+0-33 




SOI/SOAR 


9,r,i 


ACT-CL J0509-5341 


0.36 ±0.13 


n op:+0.35 


0.461 


MOSAIC/Blanco 


9,r,i 


ACT-CL J0516-5430 


0.25 ± 0.08 


26+0-^^ 


0.294 


MOSAIC/Blanco 


9,r,i 


ACT-CL J0528-5259 


0.69 ±0.05 


69+0-^'^ 


0.768 


MOSAIC/Blanco 


9ir,i 


ACT-CL J0546-5345 


0.93 ±0.13 


95+0-^^ 


1.066 


MOSAIC/Blanco 


9,r,i 



1.7 X 10^ ^ Mc7) assuming a dis tribution of massive haloes 
using the 



Tinker et al. 



([2008]) mass function. Abeh 3402 
might be in this situation, which would explain the large 
positional offset between its BCG and peak SZ decrement 
and its lack of significant X-ray emission. 

4.5. The Purity of the Sample 

The extraction methodology of SZE cluster candidates 
evolved during the 2009B season and as a result the can- 
didates observed in the first run (at the NTT) came from 
a list of targets that was later revised onc e the selection 
method was more stable (see Section |3.1[) . The revised 
sample for our second run (at SOAR) in December 2009 
contained all of the high signal detections observed on the 
NTT run. This later list, from which we targeted 49 sys- 
tems, is therefore the most uniformly observed as a func- 
tion of signal-to-noise and the one best suited to estimate 
the purity of SZE candidates. We utilized the Blanco 
archival data in the 05hr field to supplement our NTT 
and SOAR imaging, although no new clusters were found 
(beyond those previously found by SPT). By purity here 
we mean the ratio of optically-confirmed clusters to SZE 
detections. In Figure [13] we show the purity from the 
sample of 49 targeted systems as a function of signal-to- 
noise as the thick and thin solid line types (these cor- 
respond to different event binning). Since we did not 
observe every cluster candidate at signal-to-noise ratios 
below 5, for reference we also show, using dashed lines, 
the fraction of candidates greater than a given signal-to- 
noise that were actually imaged in this program. 

We achieve 100% purity for signal-to-noise ratios 
greater than 5.7 where there are 10 clusters. This drops 
to a purity value of 80% for a SNR of 5.1. We targeted 
all of the SZ candidates down to a signal-to-noise of 4.8 
and for this sample the purity is 71% (17/24). Below 
this SNR value the purity of the SZ candidates drops 
sharply. Only six clusters were optically confirmed out of 



25 candidates targeted (24% purity) over the SNR range 
4.2-4.8. However, as we were unable to target all of the 
SZ candidates in this SNR range, assigning a precise es- 
timation of the purity is not possible. Another estimate 
for the purity comes from conservatively assuming that 
none of the untargeted SZE candidates is a real system, 
in which case the overall purity would be :^ 9% in the 
signal-to-noise ratio range between 4.2 and 4.8. 

An important conclusion of this analysis is that within 
our high SNR sample (>5.7), which was fully observed, 
we report no clusters above a redshift of 1.1. As we show 
next, this distribution is fully consistent with expecta- 
tions of ACDM which are based on Gaussian fluctuations 
in the primordial matter density distribution. Deviations 
from Gaussianity can perturb the expected number of 
massive clusters at high redshift. Thus the large area 
surveyed by ACT and the well-established redshift dis- 
tribution of the sample, can be used to put important 
constraints on the level of non-Gaussianity, defined by 
the dimensi onle ss coupling parameter, /nl (see Komatsu 



Spergel|[2QQl| |Lo Verde et aLl[2008[ ). This will be the 
locus of a future study. 

4.6. Mass Threshold of the Cluster Sample 

Here we describe how we estimate the approximate 
mass limit of the observed sample. We binned the num- 
ber of clusters as a function of redshift in bins of Az = 0.2 
and fitted to the redshift distribution of massive haloes 
using the [Tinker et al.| ( |2008|) mass function and cos- 
mological parameters from [Komatsu et al. ( |2009[ ), as- 
suming that all (no) clusters above (below) the thresh- 
old were detected. The number counts were assumed 
to be Poisson-distributed in each bin and the appropri- 
ate likelihood function was calculated and used as the 
figure-of-merit function for the fits. In Figure [T4| we 
show the cumulative number distribution of optically- 
confirmed ACT SZ clusters (the histograms) per square 
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Figure 12. Spectroscopic redshift versus photometric redshift for 
the sub sample of 8 ACT clusters with known redshifts. Cir- 
cles, triangles and squares represent imaging data from the MO- 
SAIC/Blanco, EFOSC/NTT and SOI/SOAR telescopes. Error 
bars represent 68% C.L. uncertainties on the cluster photometric 
redshift. 

degree as a function of redshift. We make two assump- 
tions here: a) that the area surveyed is 455 deg^ even 
though all of this was not covered to the same depth, 
and b) that the sample is complete out to a redshift of 
1.2. The upper histogram is for the entire sample of clus- 
ters from Table 2, while the bottom histogram is for the 
subset of clusters above a SNR of 5.7 which has 100% 
purity. We overplot curves that give the expected dis- 
tribution for different mass thresholds bas ed on struc- 
ture formation in a ACDM model (Tinke r eFaL] [2008). 
The full sample of 23 is well characterized by a mass 
threshold of 8 x 10^^ Mq and the 10 highest significance 
clusters, where all ACT SZ candidate clusters were op- 
tically confirmed, are consistent with a mass threshold 
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Figure 13. The SZE sample purity greater than a given SNR 
in the ACT 148 GHz maps. The purity is defined as the ratio 
of the number of confirmed clusters to the number of observed 
clusters. The thick solid line represents the purity binned in n=2 
events, while the thin line represents no binning. The dashed line 
represent the fraction of candidates from the sample greater than 
a given SNR that were actually observed. The inset plot shows 
the cluster cumulative distribution as a function of SNR for the 
optically-confirmed cluster sample. 
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Figure 14. Cumulative distribution of ACT SZE clusters per 
square degree as a function of redshift. The top histogram is for 
the full sample of 23 clusters; the bottom one is for the high sig- 
nificance subsample of 10 clusters. The curves show the best fit 
model distributions for each. 

of 1 X 10^^ M0. These mass thresholds are in agreement 
with their X-ray luminosities (Fig. 15) and with simu- 
lations given the current ACT survey noise (Sehgal et 
al. 2010). Thus there is consistency between the SZ, op- 
tical, and X-ray observations with the model. Note that 
if the actual area were significantly smaller than 455 deg^ 
or clusters were missed, the mass thresholds would de- 
crease. Therefore the mass threshold values quoted here 
can be considered strong upper limits. 



4.7. Unconfirmed High Signal-to-Noise Candidates 

Beyond the reported sample of 23 clusters (Table [2|, 
there are in addition three SZE cluster candidates (not 
reported) at relatively high SNR (>5.1) for which we do 
not find an optical counterpart using the procedure de- 
scribed above. One of these we observed at the NTT 
and then again for twice as long at SOAR, where we 
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also included the 2:-band. Although there are potential 
BCGs in the field of view (based on our low accuracy 
photometric redshifts), none is accompanied by a red 
sequence containing more than 15 galaxies within a 
projected radius of 500 kpc. This value is consistent 
with the average surface density of galaxies ai z 1 
and additionally is far below the typical observed rich- 
ness of the confirmed clusters {c^ 60 members) in the 
optical. Thus we can confidently conclude that these 
three candidates are not optical clusters in the redshift 
range 2: < 1.2 with optical properties like those of the 
confirmed sample (Table [2]). An obvious possibility is 
that these are high redshirt clusters, beyond t he reach of 
our obser vational program. However, for the [Tinker et] 
al. (2008) distribution of massive haloes in the ACDM 
cosrnology, we expect only 0.6 or 0.1 clusters (for M = 6 
and 8 x 10^^ M©, respectively) beyond a redshift of 1.2 
over the area surveyed. Future near-IR or warm phase 
Spitzer IRAC imaging is required to determine whether 
these are high-redshift massive clusters. 

5. CLUSTER X-RAY PROPERTIES 

In order to estimate masses we turn to X-ray observa- 
tions, using archival datasets from ROSAT, Chandra and 
XMM-Newton, and measure X-ray luminosities to serve 
as coarse mass proxies. Here we describe the X-ray data 
reduction and analysis and present the X-ray properties 
of the ACT clusters. 

5.1. ROSAT Observations 

We searched for X-ray counterparts to the ACT clus- 
ters using th e RASS data followin g the same procedure 
as in Menanteau & Hughes (2009). The raw X-ray pho- 
ton event lists and exposure maps were downloaded from 
the MPE ROSAT Archiv and were queried with our 
own custom software. At the position of each cluster, 
RASS count rates in the 0.5 — 2 keV band (correspond- 
ing to PI channels 52-201) were extracted from within 
radii of 3' to 9' (chosen appropriately for each cluster) for 
the source emission and from within a surrounding annu- 
lus (12' to 30' inner and outer radii) for the background 
emission. Obvious point sources were excluded from 
both source and background regions. The background- 
subtracted count rates were converted to X-ray luminos- 
ity (in the 0.1-2.4 keV band) assuming a thermal spec- 
trum (/cT = 5 keV) and the Galactic column density 
of neutral hydrogen appropriate to the source position, 
using data from the Leiden/ Argentine/Hi Bonn survey 
Kalberl a et al | |2005). In Table H we show the X-ray 
'uxes and luminosities for all ACT clusters, regardless 
of the significance of the RASS detection. 

5.2. Chandra Observations 

All available archival Chandra data for clusters in our 
sample are used in this analysis. The combined expo- 
sure times of the observations for each cluster are shown 
m Table [5j Data are reduced using; CI AO version 4.2 
and calibration data base version 4.2.0. The data are 
processed starting with the level 1 events data, removing 
the cosmic ray afterglow correction, and generating a new 
bad pixel file that accounts for hot pixels and cosmic ray 

ftp://ftp.xray.mpe.mpg.de/rosat/archive/ 



afterglows. Using the newly generated bad pixel file, the 
charge transfer inefficiency correction, time- dependent 
gain adjustment, and other standard corrections are ap- 
plied to the data. The data are filtered for ASCA grades 
0, 2, 3, 4, 6 and status=0 events and the good time in- 
terval data provided with the observations are applied. 
Periods of high background count rate are excised using 
an iterative procedure involving creating light curves in 
background regions with 259 s bins (following the ACIS 
"Blank-Sky" Background File reduction), and excising 
time intervals that are in excess of 3 a (=rms) from the 
mean background count rate. This sigma clipping pro- 
cedure is iterated until all remaining data lie within 3 
a of the mean. All observations are in the ACIS-I con- 
figuration and background regions are chosen among the 
chips free from cluster emission, at approximately the 
same distance from the readout as the cluster. This is to 
account for a gradient in the front-illuminated chip re- 
sponse seen when invest igating deep blank-sky exposures 
( [Bonamente et al.||2004 [). 

Images are created by binning the data by a factor of 
four, resulting in a pixel size of 1.97'^ Exposure maps are 
constructed for each observation at an energy of 1 keV. 
Images and exposure maps are combined in cases where 
there are multiple observations of the same cluster. A 
wavelet-based source detector is used on the combined 
images and exposure maps to generate a list of potential 
point sources. The list is examined by eye, removing 
bogus or suspect detections, and then used as the basis 
for our point source mask. 

Spectra are extracted in regions that are 3<j above 
the background level, where a is the uncertainty in the 
background measurement. Surface brightness profiles are 
constructed using 5^' bins and are used to estimate the 
background surface brightness and the uncertainty in the 
background. The radius corresponding to 3a above the 
background, Osa^ is computed and used for the spectral 
extraction, accounting for the point source mask. Spec- 
tra are extracted and responses calculated for each indi- 
vidual observation. All of the spectra for a given cluster 
are then fit simultaneously. 

ACIS quiescent backgrounds are used to construct 
background spectra for each of the observations. Ap- 
propriate quiescent background files are reprojected to 
match each observation and background spectra are ex- 
tracted in the same Osa- regions used for the cluster spec- 
tra. The background spectra are then normalized to 
match the count rate of the cluster spectra in the 9.5- 
12.0 keV range, where the Chandra response is negligable 
an d the p article background dominates (e.g., Vikhlinin et 



aL]|2005|. 

XSFEC (|Arnaud||1996D is used to fit the X-ray spec- 
trum assuming a thermal plasma e mission model (using 
the so-called "Mekal" model; ^Mewe et al. 1985) , multi- 
plied by a Galactic absor ption m odel, witn solar abun- 
dances of Asplund et aD ('2009') and cross sections of 
Balucinska- Church & McCammon (1992) including an 
updated He cross section (Yan et al. 1998). For the ab- 
sorption component we fixed the column density of neu- 
tral hydrogen to the Galactic value toward the cluster po- 
sition (obtained in the same way as for the RASS data). 
Spectra are grouped such that there are a minimum of 20 
counts per bin and the statistic is used for the fit. All 
spectra for a given cluster are fit simultaneously to the 
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Figure 15. The soft X-ray flux (in the observed band of 0.1-2.4 keV) as a function of redshift for ACT SZE-selected clusters. The flUed 
squares are known clusters and the open circles are RASS fluxes for seven new ACT clusters with only photometric redshifts. Red circles 



and blue stars represent flux values from Chandra and XMM-Newton observations respectively. The solid and dashed curves represen t the 
expected X- rav fluxes for clusters with masses M200 = 4 x IQ-'^'^ Mq and 1 x IQ-*^^ Mq, respectively, using scaling re lations froni i Vikhlinin l 
|et al.|p009| . The gray area at the top corresponds approximately to the regime of the RASS bright source catalog ( |Voges et al.|iUU9] r 



same plasma model, including linking the normalizations 
between data sets. The fit is limited to photons within 
the energy range 0.3-9.0 keV. To estimate fluxes in the 
0.1-2.4 keV ROSAT hdiiid^ the energy range is extended 
below Chandra'^s nominal energy range. Fluxes, Fx, and 
luminosities, Lx, in the 0.1-2.4 keV band are then com- 
puted from the best fit model. Uncertainties in the fluxes 
and luminosities are estimated by varying the tempera- 
ture, Tg, abundance, Z, and normalization of the model 
using the extremities of the 68% confidence regions for 
the parameters. The results are summarized in Table [Sj 
which shows the redshift (z), total effective exposure time 
{texp)^ column density (Nh)^ the spectral extraction ra- 
dius in angular (^3^) and physical (Rsa) units, spectral 
parameters from the fit (Tg and Z), and the resulting 
flux and luminosity (Fx and Lx) for each cluster. 

5.3. XMM-Newton Observations 

For the XMM-Newton observations we started with 
the pipeline-processed data supplied as part of the stan- 
dard processing. A combined background-subtracted 
and exposure-corrected image was created for each clus- 
ter in the 0.5-2 keV band by merging data from the three 
EPIC cameras: PN, MOSl and M0S2. These were used 
to define the optimal radii for spectral extraction, ini- 
tially taken as the radius where the estimated cluster 
emission fell below the background level. These source 



regions were further refined by increasing the extraction 
radii by factors of 10% until the extracted flux values 
converged to an accuracy of a few percent. The regions, 
in nearly all cases, were elliptical. Background spectra 
came from a surrounding annular region with sufficient 
area to obtain good photon statistics. Event lists were 
filtered to remove time periods of high count rate back- 
ground flares and to include only good event patterns 
(< 12 for MOS and < 4 for the PN). Other standard 
filters to remove bad pixels and chip edges, for example, 
were applied. We used standard XMM-Newton software 
tools to calculate the instrumental response functions for 
each cluster (treating the data from each camera inde- 
pendently). The cluster image itself was used for the 
weighting to generate the response functions. The same 
absorbed thermal plasma model as described above was 
used to fit the PN, MOSl, and M0S2 data, again us- 
ing XSPEC. The individual spectral data sets were fit- 
ted separately with linked parameters. X-ray fiuxes and 
luminosities were determined using the best-fitted spec- 
trum. The uncertainties on these include the la range 
on temperature and overall normalization. Results are 
given in Table [6] 

5.4. Summary of X-ray Properties 

The soft X-ray fluxes as a function of redshift are shown 
in Figure [15] We plot all datasets: RASS in gray, Chan- 
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Table 4 

ROSAT Cluster Properties 









ti 


e 


R 






± uescriptor 


z 




^^±u cm ) 


(circmin) 


[tl^Q KpCj 


y\J.l — Z . 4t Ke V J 


{r\ -\ 9/1 l^c.A/'^ 
y\J.l — Z.4t Ke V J 


ACT-CL-J0145-5301 


0.118 


247 


2.67 


9.0 


1151 


67.59 ±8.39 


2.30 ±0.29 


ACT-CL-J0641-4949 


0.146 


859 


3.37 


6.0 


920 


6.32 ± 1.82 


0.34 ±0.10 


ACT-CL-J0645-5413 


0.167 


496 


5.60 


9.0 


1542 


91.17 ± 7.02 


6.48 ±0.50 


ACT-CL-J0638-5358 


0.222 


834 


6.06 


8.0 


1717 


69.57 ±4.74 


9.13 ±0.62 


ACT-CL-J0516-5430 


0.294 


380 


2.05 


7.0 


1844 


34.83 ±4.90 


8.43 ±1.19 


ACT-CL-J0658-5557 


0.296 


470 


4.90 


7.0 


1853 


80.26 ±6.53 


19.72 ±1.61 


ACT-CL-J0245-5302 


0.300 


621 


3.15 


7.0 


1870 


23.39 ± 3.54 


5.92 ±0.90 


ACT-CL-J0217-5245 


0.343 


357 


2.68 


7.0 


2047 


30.56 ±4.86 


10.40 ±1.65 


ACT-CL-J0237-4939 


0.40 


778 


2.55 


3.5 


1128 


3.84 ± 1.28 


1.84 ±0.61 


ACT-CL-J0235-5121 


0.43 


460 


3.03 


3.5 


1178 


10.91 ± 2.44 


6.14 ± 1.37 


ACT-CL-J0707-5522 


0.43 


620 


6.37 


3.5 


1178 


11.29 ±2.28 


6.36 ±1.28 


ACT-CL-J0330-5227 


0.440 


746 


1.44 


6.0 


2046 


35.77 ±3.35 


21.20 ±1.98 


ALy 1 -LyL/-J UoUy-Oo4i 


0.461 


418 


1.46 


Q C 
O.O 


1224 


A AO _L 1 

4.4c) ± 1.7b 


o on _L 1 1 
z.yU ± 1.15 


ACT-CL-J0304-4921 


0.47 


699 


1.77 


5.0 


1770 


8.15 ± 1.89 


5.61 ± 1.30 


ACT-CL-J0215-5212 


0.51 


369 


2.49 


5.0 


1850 


3.79 ±2.28 


3.13 ±1.88 


ACT-CL-J0438-5419 


0.54 


308 


1.02 


3.5 


1334 


21.48 ±3.83 


20.21 ±3.60 


ACT-CL-J0346-5438 


0.55 


774 


1.44 


3.0 


1154 


3.58 ± 1.13 


3.52 ±1.11 


ACT-CL-J0232-5257 


0.59 


425 


2.84 


3.5 


1393 


2.94 ± 1.59 


3.39 ±1.83 


ACT-CL-J0559-5249 


0.611 


995 


5.06 


4.0 


1618 


4.69 ± 1.27 


5.85 ±1.58 


ACT-CL-J0616-5227 


0.71 


1223 


4.24 


4.0 


1724 


5.08 ± 1.28 


8.92 ±2.24 


ACT-CL-J0102-4915 


0.75 


227 


1.72 


4.0 


1761 


9.25 ±3.27 


18.39 ±6.49 


ACT-CL-J0528-5259 


0.768 


718 


3.24 


3.0 


1332 


1.26 ±0.95 


2.64 ±2.00 


ACT-CL-J0546-5345 


1.066 


1092 


6.78 


3.0 


1462 


2.04 ±0.88 


9.14 ±3.94 


^ Units are xlO"-'^'^ erg 


— 1 —2 

s cm . 















^ Units are x 10^^ erg s ^ . 

dra in red, and XMM-Newton in blue. We indicate the 
region (at the high-end flux) that approximately corre- 
sponds to the flux h mit of the RASS bright source catalog 
(Voges et al. 1999). In cases of multiple X-ray datasets 
the agreement among the observed fluxes is good, with 
the possible exception of ACT-CL-J0330-5227, where the 
emission from the foreground cluster contaminates the 
RASS flux estimate. 

For reference we show curves of the expected 
(observed-frame) X-ray fluxes for clusters with assumed 
masses of M200 = 4 x 10^"^ M© (dashed) and 1 x lO^^M© 
(solid) using the X-ray lum inosity versus mass scaling re- 
lation in Vikhlinin et al.| (|2QQ9). We use their eqn. 22, 
which includes an empirically determined redshift evolu- 
tion. We convert their X-ray band (emitted: 0.5-2 keV 
band) to ours (observed: 0.1-2.4 keV) assuming a ther- 
mal spectrum at the estimated cluster temperature de- 
termined usings; the m ass-temperature relation also from 
lyikhlinin et al. (|2009). This too has a redshift depen- 
dence, so the estimated temperatures vary in the ranges 
2.9-5.0 keV and 5.2-9.0 keV over 0.0 < z < 1.1 for the 
two mass values we plot. We use conversion factors as- 
suming the redshift-averaged temperatures of 4 keV and 
7 keV, since the difference in conversion factor over the 
temperature ranges is only a few percent, negligib le on 
the scale of Figure [15] The mass values in Vikhlinin et| 



to 1.7 over 0.0 < z < 1.1, and depends only weakly on 
cluster mass (few percent) for the 2 values plotted here. 

The X-ray fluxes of the ACT SZE-selected clusters 
scatter about the 1 x 10^^ Mq curve, validating the mass 
threshold we estimated above from the cluster number 
counts. The most extreme outlier (on the low flux side) 
is Abell 3402, indicating a considerably lower mass for 
this system, assuming (as we have) that the SZE candi- 
date is at the redshift of Abell 3402. For those clusters 
with archival Chandra or XMM-Newton data our mea- 
sured temperatures are typically high kT > 5 keV and 
in good agreement with results f rom the rec ent SPT clus- 
ter X-ray analysis (Andersson et al. 2010). The masses 



inferred from the X-ray temperatures are consistent with 
the mass curves shown in Figure 15 



al 



(2009) are deflned with respect to an overdensity of 
times the critical density of the Universe at the clus- 
ter redshift, while in this paper we deflne cluster masses 
with respect to an overdensity of 200 times the average 
density at the cluster redshift. We calculate the conver- 
sion factor assuming an NFW halo mass proflle with the 
power-law r elation between con centration and mass de- 
termined bylDuffy et al. (2008) from large cosmological 
A^-body simulations. I'his mass conversion factor is ap- 
proximately 1.8 averaged over redshift, varies from 2.2 



6. CONCLUSIONS 

After many years of intense effort, we have presented 
the flrst optically-conflrmed sample of 23 SZE-selected 
clusters of galaxies from ACT over a cosmologically sig- 
niflcant volume. Using a combination of archival and 
targeted multi-color observations using 4-m telescopes we 
have conflrmed the presence of a BCG and a rich red se- 
quence of galaxies in 14 new SZE clusters (4 overlapping 
with previous ACT and SPT detections), of which 10 are 
presented here for the flrst time. Using a careful selec- 
tion method of SZE candidates combined with the results 
from our observations we have determined the purity of 
the ACT SZE sample to be 100% for SNR > 5.7, 80% 
for SNR > 5.1, and 71% for SNR > 4.8. For lower SNR 
values we did not target all SZE candidates; still it is 
clear that the purity of the lower signiflcance candidates 
has fallen substantially to between :^9% and 24%. 

By comparing the number distribution of clusters of 
galaxies as a function of redshift with the expectations 
from ACDM simulation, we estimate a typical lower mass 
limit of 8.0 x 10^"^ Mq for our fuh sample of 23. Our 
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Table 5 

Chandra Cluster Properties 







texp A^iJ Oscr R3(T 


Te 


Z 


Fx"" Lx"" 


ACT Descriptor 


z 


(ks) (10^0 cm-2) (arcmin) (/i~^ kpc) 


(keV) 




(0.1-2.4 keV) (0.1-2.4 keV) 



ACT-CL J0638-5358 


0.222 


19.9 


6.06 


5.20 


1116 


ACT-CL J0516-5430 


0.294 


9.5 


2.05 


4.05 


1068 


ACT-CL J0658-5557 


0.296 


589.4 


4.90 


5.92 


1568 


ACT-CL J0330-5227 


0.440 


19.3 


1.44 


3.60 


1229 


ACT-CL J0509-5341 


0.460 


28.8 


1.46 


1.95 


683 


ACT-CL J0528-5259 


0.768 


27.2 


3.24 


1.04 


463 


ACT-CL J0546-5345 


1.068 


57.1 


6.78 


1.24 


605 



8.51^ 



+ 0.36 
0.34 

9.59ii:g^3 

13.56ig:ll 
4_32+o.2i 



-0.19 
□ +1.43 
-1.07 

4,67+1.45 



9.39^ 



8.54' 



0.99 
+ 1.38 
1.05 



46+°-°^ 

U.^D_Q Qg 

0,49+0.25 



-0.24 

r. + 0.02 

" -0.02 
Q_-Ll+0.07 



0.36^ 



0.35 



0.33 



0.07 
+ 0.22 
0.23 

0.30 

+ 0.29 
0.29 



76.64 
26.79 
87.75 
25.72 
6.61 
1.02 
2.06 



1.02 
1.02 
1.00 
0.93 
+0.22 
0.22 
+0.62 
0.60 
+0.26 
0.24 
+0.12 
0.10 
+0.12 
0.11 



9.87^ 



7+0.13 
-0.13 

6.30lgg 
20.69ig;0^ 
15.50^°-3" 



-0.36 
^+0.16 
-0.15 
2,15+0.24 



4.14^ 



8.43 



0.22 
+0.47 
0.44 



^ Units are xlO erg s 
^ Units are xlO^'^ erg s-^ 



Table 6 

XMM-Newton Cluster Properties 



ACT Descriptor 


z 


(ks) 


Nh 
(1020 cm-2) 


^minor 5^major 

(arcmin) 


D b 
J^mean 

{h^o kpc) 


Te 

(keV) 


Z 


Fx"" 
(0.1-2.4 keV) 


Lx^ 
(0.1-2.4 keV) 


ACT-CL J0145-5301 


0.118 


37,22 


2.67 


6.5,10.8 


1031 


5.60 ±0.08 


0.30 ±0.02 


82.96 ±0.52 


2.83 ±0.02 


ACT-CL J0645-5413 


0.167 


44,28 


5.60 


5.4,8.9 


1188 


7.49 ±0.09 


0.24 ±0.01 


118.70 ±0.40 


8.36 ±0.03 


ACT-CL J0516-5430 


0.294 


11,8 


2.05 


6.2,8.0 


1856 


7.44 ± 0.38 


0.21 ±0.05 


37.57 ±0.47 


8.93 ±0.11 


ACT-CL J0658-5557 


0.296 


33,21 


4.90 


3.9,5.1 


1181 


10.80 ±0.22 


0.22 ±0.03 


65.22 ±0.39 


15.51 ±0.09 


ACT-CL J0330-5227 


0.440 


71,57 


1.44 


4.0,4.0 


1365 


5.46 ±0.27 


0.30 


24.28 ± 2.34 


14.60 ± 1.40 


ACT-CL J0559-5249 


0.611 


20,14 


5.06 


2.0,3.5 


1070 


8.09 ±0.75 


0.30 


5.10 ±0.10 


6.06 ±0.12 



^ The quoted MOS exposure time is the average of the two MOS cameras. 
^ Geometric mean of the major and minor axes. 
^ Units are xlO"-*^^ erg s"-*^ cm~2. 
^ Units are xlO^^ erg s-^ 



highest significance subsample of 10 ACT SZE cluster 
candidates (above a SNR of 5.7) were all confirmed opti- 
cally is consistent with a mass threshold of 1.0 x 10^^ Mq. 
We present the results of analysis of archival X-ray data 
from ROSAT, Chandra, and XMM-Newton. The X-ray 
data yield uniformly high temperatures {kT 2^ 7 — 17 
keV) and luminosities {Lx — 10^^ — 10^^ erg s~^), which 
are broadly consistent with the mass thresholds inferred 
from the cluster number counts as a function of redshift. 
Further work on measuring cluster masses using X-ray 
data and other techniques (e.g., galaxy dynamics and 
weak lensing) is underway. 

While the ACDM model predicts that there should be 
very few clusters with masses greater than 8 x 10^^ Mq 
at 2; > 1.2, (i.e., 0.1 clusters with our current area and 
ma ss limit or 0.4% of the sample) early dark en ergy mod- 
els ( Bartelmann et ar]|2006 Alam et al.||2010 ) an d mod- 



els with sig niticant primordial non-Gaussianity (Cayon 
et al.|[2QTQ ^) predict significantly more high redshitt clus 



ters. Moreover, in quintessence models, the number of 
high redshift clusters is also a sensitive function of the 
dark energy equation of state ( Alimi et al.|2010 |). There- 
fore, the detection of even a single massive high redshift 
cluster would challenge our current standard paradigm. 
However our ability to investigate this question here is 
constrained somewhat by the upper redshift limit set by 
the depth of our optical program. Future follow-up IR 
observations of high significance SZE candidates that are 
currently undetected in our optical imaging can probe 
the existence of z > 1.2 massive systems and provide a 



simple yet powerful test of ACDM. Even in the event of a 
null result (i.e., none of the high significance SZE decre- 
ments turns out to contain a cluster), the ACT cluster 
sample has the potential to set strong constraints on al- 
ternative dark energy models. 

The observations on which this paper were based rep- 
resent the marriage of two different communities (CMB 
and optical) in multiple countries working for a com- 
mon goal. In particular, the optical observations were 
coordinated and led by Felipe Barrientos and Leopoldo 
Infante (Pontificia Universidad Catolica de Chile) and 
John P. Hughes and Felipe Menanteau (Rutgers Uni- 
versity). This work was supported by the U.S. Na- 
tional Science Foundation through awards AST-0408698 
for the ACT project, and PHY-0355328, AST-0707731 
and PIRE-0507768 (award number OISE-0530095). The 
PIRE program made possible exchanges between Chile, 
South Africa, Spain and the US that enabled this re- 
search program. Funding was also provided by Prince- 
ton University and the University of Pennsylvania. We 
also acknowledge support from NASA/XMM grants 
NNX08AX55G and NNX08AX72G to Rutgers Univer- 
sity. Computations were performed on the GPC super- 
computer at the SciNet HPC Consortium. SciNet is 
funded by: the Canada Foundation for Innovation un- 
der the auspices of Compute Canada; the Government 
of Ontario; Ontario Research Fund Research Excel- 
lence; and the University of Toronto. This research is 
partially funded by "Centro de Astrofisica FONDAP" 
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15010003, Centre BASAL-CATA and by FONDECYT 
under proyecto 1085286. The observers (FM, JPH, JG, 
LI) would like to thank the La Silla, CTIO and SOAR 
staff for their support during the runs. The SOAR 
Telescope is a joint project of: Conselho Nacional de 
Pesquisas Cientificas e Tecnologicas CNPq-Brazil, The 
University of North Carolina at Chapel Hill, Michigan 
State University, and the National Optical Astronomy 
Observatory. 
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